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ABSTRACT 



We perform axisymmetric resistive MHD calculations that demonstrate that centrifugal disks can indeed form around Class objects 
despite magnetic braking. We follow the evolution of a prestellar core all the way to near-stellar densities and stellar radii. Under 
flux-freezing, the core is braked and disk formation is inhibited, while Ohmic dissipation renders magnetic braking ineffective within 
the first core. In agreement with observations that do not show evidence for large disks around Class objects, the resultant disk 
forms in close proximity to the second core and has a radius of only as 10 R a early on. Disk formation does not require enhanced 
resistivity. We speculate that the disks can grow to the sizes observed around Class II stars over time under the influence of both 
Ohmic dissipation and ambipolar diffusion, as well as internal angular momentum redistribution. 
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1. Introduction 

Understanding how protostellar and protoplanetary disks form is 
of fundamental importance to theories of star- and planet forma- 
tion. Observations show th e ir ubi quity around Class II objects 
(e.g., lAndrews & Williamsl I2005I) . In recent years, doubt was 
cast on their accepted formation mechanism, when it was shown 
that for flux freezing magnetic braking is so effective in remov- 
ing angular momentum from the parent cor e that large-scale 
(sa 10 2 AU) disks are s uppressed entirely (lAllen. Li. & Shul 
I20031 IMellon & Lil 120081 iHennebelle & Fromand. I2008I) . This 
scenario held true even when a s implified version of ambipo- 
lar diffusion (IMellon & LiL [2009) was included in the model, 
and has b een referred to as the magn etic braking catastrophe. 
Recently, Henneb elle & Ciardl (|2009) demonstrated that incli- 
nation effects can modify the efficiency of magnetic braking, 
but a supercritical mass-to-flux ratio by a factor > 3 - 5 (i.e., 
a we ak magnetic field) was still required to form a large-scale 
disk. Pu ffin & Pudritj (l2009) performed three-dimensional sim- 
ulations with ambipolar diffusion, but only resolved the first 
core, and did not find Keplerian motion. 

Runaway collapse of a prestellar core can effectively trap the 
magn etic flux in the prestellar phase (e.g jBasu & Mouschoviasl 
1994). If the evolution continued to proceed under flux-freezing, 
a big magnetic flux problem would remain, since the emerging 
star would hold 10 3 - 10 5 times more magnetic flux than ob- 
served in T Tauri stars. At densities < 10 12 crrT 3 , ambipolar dif- 
fusion causes flux leakage, while at even higher densities, matter 
decouples entirely from the magnetic field, and O hmic dissipa- 
tion becomes dominant (e.g., Nak ano et all 120021). Both effects 
are revitalized after the formation of a central star dLi & McKe4 
Il996t IContopoulos etail 1 1998b . Recently, Krasnopolskv et al.1 
(2010) have shown that for an isothermal core without self- 
gravity, only an 'anomalous' resistivity — a factor of 100 larger 
than the canonical level — allows disks of size 10 2 AU to form 
during the Class phase. However, their simulations are domi- 



nated by numerical reconnection events that make precise state- 
ments about the efficacy of magnetic braking difficult. 

Currently, there is no evidence for the presence of centrifu- 
gal d isks larger than » 5 AU around Class or Class I objects 
fe.g.. lMaurvetall l2010). However, there are outflows observed 
even at these early ages. It is therefore reasonable to assume that 
disks form at a small scale and only subsequently grow to the 
larger sizes observed in the Class II phase. We demonstrate the 
first part explicitly by using a canonical level of Ohmic dissipa- 
tion alone, and speculate that the combined effects of ambipolar 
diffusion and Ohmic dissipation will allow for the second part. 
Additionally, an initially small disk could expand significantly if 
angul ar momentum transport is regulated b y internal processes 
(e.g.. lBasuiri998tlVorobvov & Basull2007|) . 

Mach ida et al.l d2007l) performed three-dimensional simula- 
tions of resistive MHD on a nested grid, following the evolution 
to stellar densities, but were only able to integrate until a few 
days after stellar core formation. We extend their calculations in 
a dimensionally-simplified model in order to simultaneously ad- 
dress the magnetic flux problem, integrate further in time, and 
study the formation of a centrifugal disk. We show that catas- 
trophic magnetic braking can be avoided, and that a small disk 
forms in a very early phase of evolution. 



2. Method 

We solve the normalized M HD equations in axis y mmet- 
ric thin-disk geome tr y (se e ICiolek & Mouschoviasl [l993; 
Basu & Mouschovias, 1994), assuming vertical hydrostatic 
equilibrium in a vertical one-zone approximation. An integral 
method for calculating the self-gravity of an infinitesimal ly-thin 
disk is used (detailed in ICiolek & Mouschoviasl 1 19931) . with 
modifications for the finite extent and finite thickness of the flat- 
tened core. 
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In our model, the magnetic field points solely in the verti- 
cal direction inside the disk, but also possesses radial and az- 
imuthal components (B r and B^) at the disk surfaces and above. 
B r is determined from a potential field assuming force-free and 
current-free conditions in the external medium. We calculate 

and implement magnetic braking using a steady-state ap- 
proximation to the transport of Alfven w aves in the external 
medium, as in lBasu & Mouschoviasld 19941) . Owing to numerical 
complexity, a calibration of this method with results of three- 
dimensional MHD wave propagation through a stratified com- 
pressible medium has not been done to date. We modify the 
ideal-MHD induction equation to include Ohmic dissipation: 

dB zm ^ 1 d t \ _t] d ( dB ZiBq \ dr] dB zm 
dt rdr^ " eq ' r dr\ dr J dr dr 

Here, B zm denotes the z-component of the magnetic field at the 
midplane of the disk, and v r is the radial component of the neu- 
tral velocity. 

We use the paramet rization of iMachida et al.l d2007l) for the 
resistivity calculated bv lNakano et al.l {2002), with a dimension- 
less scaling parameter rjo whose standard value is unity. The re- 
sistivity is then 

x l-tanh( — -| cm 2 s _1 , (2) 

\10 15 cm- 3 /J V ' 

where n is the volume number density, and the term in square 
brackets is a cutoff representing the restoration of flux freez- 
ing at high densities. Th e uncertainties in rjp h inge largely on 
the grain prop e rties ( e.g. jMachidaet all 12007). Different from 
Mac hida et al.l d2007f) . we do not (inconsistently) pull the resis- 
tivity outside all spatial derivatives. 

For simplicity, we replace the detailed energy equation 
by a baro tropic relati o n. Th e temperature-density relation of 
Masunaga & Inutsukal ([2000) is transformed into a pressure- 
density relation using the ideal gas law P = nk^I ', where P is the 
pressure, k% is Boltzmann's constant, and T is the temperature. 
We calculate the midplane pressure self-consistently, including 
the effects of the weight of the gas column, constant external 
pressure (P e xt = 0.1 jiGI^/I), magnetic pressure, and the extra 
squeezing added by a central star (once present). 

T he MHD equat ions are solved with the method of lines 
(e.g., ISchiesserl 1 19911) using a finite volume approach on an 
adaptive grid with up to 1024 radial cells in logarithmic spacing. 
The smallest cell is initially 10~ 2 AU and as small as 0.02 R G at 
the highest refinem ent. We use the second-order van-Leer TVD 
advection scheme dvan Leerlll977l) . and calculate all derivatives 
to second-order accuracy on the nonuniform grid. The code will 
be described in detail in a forthcoming paper. 

3. Initial conditions and normalization 

We assume that our initial state was reached by core 
contraction preferentia ll y alo ng magnetic field lines (e.g., 
iFiedler & Mouscho viasl 1 1993b and rotational flattening, and 
start with initial profiles for the column density and angular ve- 
locity given by 

So 2Q C 
2(r) = -=?==, Q(r)= — ■ (3) 

VI + (r/R) 2 VI + (r/R) + 1 



Here, R « 1 , 500 AU approximately equals the Jeans length at 
the core's initial central density (see below). The column den- 
sity profile i s representative of th e early stage of collapse (e.g., 
iBasul 119971: iDapp & Basul 120091) . and the angular velocity pro- 
file reflects that the specific angular momentum of any parcel is 
proportional to the enclosed mass. 

We assume an initial profile for B zeq in a way that the nor- 
malized mass-to-flux ratio jj. = S/B zeq 2k VG = 2 everywhere, 
whic h is the approximate s tarting point of runaway collapse 
(e.g., iBasu & Mouschoviasl 1 19941) . The radial velocity is ini- 
tially zero. The initial state is not far from equilibrium, because 
the pressure gradient and magnetic and centrifugal forces add up 
to w 82% of the gravitational force. Our results do not depend 
strongly on the choice of initial state as long as gravity remains 
dominant. 

The initial central column density and number density are 
So = 0.23 g crrT 2 and n c = 4.4 x 10 6 cm 4 , respectively. The 
total mass and radius of the core are 2.5 M and 1.2 x 10 4 AU, 
respectively. The initial central magnetic field strength is B zfiq * 
200 fiG. We choose the external density in a way that n c /n ext = 
500, (i.e., n ext =3 10 cm ), and the central angular velocity Q c 
so that the cloud's edge rotates at a rate of 1 km s " 1 pc" 1 , consis- 
tent with observations of molecular cloud cores ( Goodman et al.L 
ll993llCaselri et al.Ll2002l) . 



4. Results 

4. 1. Prestellar phase and formation of the second core 

During the prestellar phase (for number densities n < 10 11 cm -3 ) 
the collapse proceeds in a nearly self-similar fashion. We find 
that — insensitive to initial conditions — the column density is 
approximately oc r~ x for three orders of magnitude of central 
enhancement, which corresponds to the volume density being 
oc r~ 2 for a central enhancement of a; 10 6 . This profile i s char - 
acteristic of a collapsing prestellar core (e.g., iLarsonl 119691) . 
The collapse proceeds dynamically, and to a good approxima- 
tion under isother mality, flux-freezing, and w ithout significant 
magnetic braking dBasu & Mousc hoviasl 1 1 994b . 

Once the density reaches « 10 11 cm~ 3 , the central region 
becomes opaque and traps the energy released by the collapse, 
which pre viously could e scape freely as radiation. Th is region 
heats up dLarsonl 119691: iMasunaga & Inutsukal 120001) and its 
thermal pressure gradient temporarily stabilizes it against fur- 
ther collapse. This is theirs? core. Its density and temperature 
increase with continued accretion, while its size stays almost 
constant at a; a few AU, bounded by an accretion shock. The 
external gravitational potential of this object closely resembles 
that of a point mass, and an expansi on wave dev elops and moves 
outward at nearly the sound speed dShulll977l) . Material within 
this region moves at near free-fall speed. 

When the temperature in the first core reaches * 2000 K, for 
n > 10 15 cm~ 3 , hydrogen molecules are collisionally dissoci- 
ated. This process provides an energy sink, so that the tempera- 
ture rise stagnates, and the collapse reinitiates. As the tempera- 
ture rises yet further, hydrogen is ionized sufficiently that flux 
freezing is re-established. Collapse is then finally halted, and 
sufficiently high de nsities are reached that electr on degeneracy 
becomes important (Masunaga & Inutsuka, 2000). A protostel- 
lar core (the s econd core) forms with a radius w a few R (e.g., 
ILarsonl [l969). This Class object initially only has a mass of 
a few xl0~ 3 M . The gravitational potential resembles that of 
a point mass outside the second core, and an expansion wave 
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once again moves outward from the accretion shock, eventually 
consuming the entire region of the previous first core. 

Figure Q] shows the profiles of column density, mass-to-fiux 
ratio and angular velocity shortly after the second core forms (« 
4.8 x 10 4 yr into the simulation). For n > 1 12 cm' 3 , Ohmic dis - 
sipation becomes dynamically important (Nakano et al., 2002), 
because all charge carriers decouple from the magnetic field, and 
flux is dissipated. While the density in the first core increases, 
we find th e magnetic field strength remains stagnant. A mag- 
netic wall dLi & McKeeL [1996L IContoPoulos et all fl998h forms 
at » 10 AU, visible as a sharp transition in column density in 
the resistive model (rjo = 1, top panel). Here, infalling neu- 
trals within the expansion wave are temporarily slowed down by 
the relatively well-coupled magnetic field that is expelled from 
the first core with a radius * 1 AU. Further inward, the neu- 
trals resume near-free-fall motion, but with enhanced magnetic 
support and at a greater column density than for flux-freezing 
(770 = 0, dotted line). Under angular momentum conservation 
(no magnetic braking), the additional rotational support stabi- 
lizes the first core against further collapse (to p panel, dash-dotted 
line), consistent with previous findings (e.g. JSaigo & Tomisakal 

1200a) . 

Because of magnetic flux dissipation, the mass-to-flux ratio 
increases by almost three orders of magnitude in the first core re- 
gion for 770 = 1, but by almost two orders of magnitude even for 
770 as low as 0.01 (Fig.Q] middle panel). The torque on the cloud 
caused by mag netic braking scales linearly with the amount of 
enclosed flux (Basu & Mous choviasl fl994). Ohmic dissipation 
therefore allows spin-up to proceed, even though the rotation 
rate is still reduced by a factor of a few outside the first core, 
compared with the case without magnetic braking (Fig. Q] bot- 
tom panel, dash-dotted line). In the flux-freezing case, the com- 
paratively slow evolution of the first core allows enough time for 
magnetic braking to spin down the first core region, and 'catas- 
trophically' brake it (Fig.Q] bottom panel, dotted line). 

4.2. Evolution after second core formation 

When the second core forms, the thin-disk formulation breaks 
down, because the object is now truly hydrostatic and spherical. 
Presumably, dynamo processes within the fully convective pro- 
tostar will also take over, and the magnetic field will mostly de - 
couple from that of its parent core (Me stel & LandstreetL 12005). 
Therefore, we switch off magnetic braking in the second core, 
and introduce a sink cell with a size of 3 R e , slightly larger than 
the second core. The processes within it are beyond the scope 
of our model, but are not expected to significantly influence 
the surroundings. This is not necessarily the case with a sink 
cell of size » 10 AU, as is the more common ap proach (e.g., 
IVorobvov & Basull2007tlMellon & LiLl2008l 12009b . 

Figure|2]shows the profiles of column density, infall velocity, 
and the ratio of centrifugal to gravitational acceleration about a 
year after the introduction of the sink cell. Centrifugal balance 
is achieved in a small region (« 10 R e ) close to the center (bot- 
tom panel) in the resistive model. This is a necessary and suf- 
ficient condition for the formation of a centrifugally-supported 
disk. At the same time all infall is halted there and the radial ve- 
locity plummets (middle panel). After a few years of evolution, 
a Toomre instability develops, and the rotationally-supported 
structure breaks up into a ring (top panel). At this point, we stop 
the simulation, because more physics would be required to fol- 
low the further evolution of the disk. Our model allows a clear 
distinction between a magnetic pseudo-disk, a flattened (disk- 
like) prestellar core, and a centrifugal (nearly Keplerian) disk. 



10~ 4 10~ 2 10° 10 2 10 4 




10" 4 10~ 2 10° 10 2 10 4 

r / AU 

Fig. 1. Spatial profiles of various quantities after the second col- 
lapse (after m 4.8 x 10 4 yr). Top: The first and second core and 
their accretion shocks are at radii m 1 AU and « 5 X 10~ 3 AU ~ 
1 R Q , respectively. Within the expansion wave outside the first 
core, the column density profile assumes that of free-fall col- 
lapse in the flux-freezing case (770 = 0), and shows a magnetic 
wall in the resistive case. Beyond w 20 AU, the prestellar in- 
fall profile remains unchanged. Without magnetic braking (dash- 
dotted line), the first core is larger and rotation prevents further 
collapse. Middle: The mass-to-flux ratio is increased by (even 
weak) Ohmic dissipation by > 10 2 . The influence is significant 
even well outside the boundary of the first core (at a few AU). 
Bottom: For flux-freezing, catastrophic magnetic braking spins 
down the first core to nearly the background rotation rate. In the 
resistive case (solid line), the rotation rate outside the first core is 
reduced only slightly compared with the case without magnetic 
braking (dash-dotted line). 



This distinctio n is not clear in profiles from three-dimen sional 
simulations dMachida et all l2007t iDuffin & Pudritzll2009h . 

Figure [3] shows the magnetic field line topology above and 
below the disk on two scales (10 AU and 100 AU), for both flux- 
freezing and resistive models. They are calculated immediately 
after the formation of the second core, as suming force-free an d 
current-free conditions above a thin disk (Meste l & Ravi [19851) . 
The split monopole of the 770 = model (dashed lines) is cre- 
ated as field lines are dragged in by the freely falling material 
within the expansion wave front at * 20 AU. This is replaced 
by a much more relaxed field line structure in the resistive case 
(solid lines). The extreme flaring of field lines in the 770 = 
mode l is a fundamenta l cause of the magnetic braking catastro- 
phe. |GalUetal] d2009l) presented similar field configurations re- 
sulting from a simplified model for resistive collapse. 
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Fig. 2. Spatial profiles of various quantities * 1 yr after the 
introduction of a sink cell of size a 3 R G . Top: The Toomre- 
unstable centrifugally-supported disk breaks up into a ring. 
Middle: Infall is halted after the formation of a centrifugal disk 
at « 5 X 10~ 2 AU ~ 10 R e in the resistive case (rjo = 1), while for 
flux-freezing (rjo = 0), infall continues. Bottom: Ratio between 
centrifugal and gravitational accelerations. The dashed line in- 
dicates rotational balance, achieved within » 10 R e with Ohmic 
dissipation. For flux-freezing, rotational support is negligible in 
the first core region owing to the magnetic braking catastrophe. 





Fig. 3. Magnetic field lines. The box on the left has dimensions 
10 AU on each side, while the box on the right has dimensions 
100 AU. The dashed lines represent the flux-freezing model 
(fjo = 0), while the solid lines show the same field lines for the 
resistive model (rjo = 1). The second core has just formed and is 
on the left axis midplane. 



probe for disks down to « 10 AU. We anticipate that the centrifu- 
gal disk that forms in our simulations can grow over time into 
disks of size » 100 AU ob served around Class II objects. Recent 
work (Machid a et all I20T0T) shows that magnetic braking can be 
cut off at late times as the envelope is accreted, and the existing 
disk can also grow by internal angular mo mentum redistribu- 
tion processes (e.g.. lVorobvov & Basull2007h. Furthermore, we 
specu late that ambipolar diffusion (e.g.. iKunz & Mouschoviasl 
1201 0b has the potential to dissipate enough flux outside the first 
core (an area not significantly affected by Ohmic dissipation) to 
reduce braking and to allow the disk to form there as well. We 
will present results of a study including both non-ideal MHD 
effects and grain physics in an upcoming paper. 
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5. Discussion and conclusions 

We demonstrate the formation of a centrifugally-supported disk 
despite the presence of magnetic braking. The magnetic braking 
catastrophe is averted by including the canonical level of Ohmic 
dissipation, which removes large amounts of magnetic flux from 
the high-density region of the first core. In the absence of Ohmic 
dissipation, this region would be spun down tremendously prior 
to the second collapse. We emphasize that disk formation hap- 
pens very shortly after the second collapse in a region very close 
to the central object, while it is still very small (< 10~ 2 M ). 
This is consistent with the observational evidence of outflows at 
a very young age. 

Our simulations yield « 0.1 — 1 kG ma gnetic fields, compa - 
rable to those observed in T Tauri stars (e.g.. |Johns-KrullLl2007h . 
in a central object of mass w 10~ 2 M Q . This is achieved by non- 
ideal MHD effects reducing the field strength by w 10 3 compared 
to a flux-freezing model. Our model does not have the capabil- 
ity of including outflows or jets, even though those are launched 
very close to the stellar surface. 

There is presently no evide nce for centr i fugal disks > 50 AU 
aroun d Class objects (e.g., lAndre et all 120021: iMaurv etail, 
l2010l) . ALMA will allow observers to improve on this, and to 
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